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ABSTRACT 


The multi-phonon processes in incoherent scattering of slow neutrons by crystals are dis- 
cussed, assuming the harmonic approximation for the erystal vibrations. The differential scat- 
tering cross section is expanded in the Hermite orthogonal functions and approximate express- 
ions for the cross section are derived. Extensive numerical calculations have been carried out 
to illustrate the accuracy of the approximations made. An approximation for the total cross 
section (the mass-ratio expansion) suggested by Placzek is discussed and in some respects 
generalized. The approximations for the differential cross section mentioned above are also 
used to derive approximate formulae for the total cross section valid for cold neutrons but 


arbitrary temperatures and mass ratios. 


Introduction 


The basic ideas of the theory of slow neutron scattering by crystals were 
developed by Wick [1], Pomeranchuk [2], Seeger and Teller [3] and Akhieser 
and Pomeranchuk [4]. A quantitative account was given by Weinstock [5], who 
discussed the temperature dependence of the total scattering. Afterwards the 
formal treatment was completed especially by Froéman [6]. He separated the 
scattering into phonon processes and consistently used the analogies with X-ray 
diffraction. An alternative method, very convenient for calculating the total 
scattering cross section, was later suggested by Placzek [7]. Recently the theory 
was reformulated by Glauber [8] and Van Hove [9] making it more surveyable. 
They derived closed expressions for the differential scattering cross section, 
which seem to be a convenient starting point for quantitative discussions. Van 
Hove also generalized the theory to general systems of nuclei, as for instance 
liquids and magnetic materials. A large number of experiments have been per- 
formed and these mainly confirm the basis of the theoretical treatment [10-19]. 
Though the formal treatment seems to be fairly complete no complete quantit- 
ative investigations have been made. Some important aspects have been thor- 
oughly discussed, as for instance one-phonon scattering [20-25] and total scat- 
tering [7, 26-29], but no satisfactory treatment of the multi-phonon processes 
is given [30, 34]. Especially in recent years great progress has been made in 
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getting detailed information on the vibrational properties of a crystal [14-18, 
31-33] and therefore it would be of importance to conclude from theory how 
the scattering should depend on the detailed properties of the crystal vibrations. 
In cases where the details give no significant effects it is important to know 
more precisely what average property of the vibrations is found from the scat- 
tering data. In this paper and an earlier one [34] these aspects of the problem 
have been discussed and where justified simplifying approximations have been 
made to get more simple expressions for the scattering cross sections. 

In part I a brief account of the theoretical background is given. In part II 
the energy distribution of the scattered neutrons in different scattering directions 
is investigated and in part III the total cross section is treated. When dealing 
with the high-phonon processes it seems to be convenient to separate the scat- 
tering into an incoherent part and an interference part, where these two should 
be treated in rather different ways. In parts II and III only the incoherent 
part of the scattering is treated. 


Part I 


GENERAL THEORY 


Chapter 1. Crystal dynamics 


A rigorous theory of the thermal vibrations in a real crystal would be very 
involved because of several physical effects which complicate the mathematical 
treatment while having little bearing on the case considered here. In order to 
obtain a simplified description of the vibrational states in the crystal we will 
therefore from the start disregard some aspects and concentrate on the im- 
portant features. We consider the crystal as perfect, apart from thermal dis- 
placements, and neglect all surface distortions of the vibrations. In addition 
we adopt the harmonic approximation for the vibrational potential. On these 
assumptions the general state of motion of the nuclei can be represented as a 
superposition of plane waves, each of which is characterized by its wave vector 
q, frequency and polarization vector C*. 

An account of elements of the dynamical theory of crystals can be found 
in the textbooks and review articles referred to [385-38]. Here we only explain 
the notations to be used and write the formulae needed in the following. As 
usual we specify the position vector of a nucleus by the lattice vector R’, 
where J=(1,,1,, 13) are the cell indices, by the equilibrium position vector in 
an elementary cell r, (k=1,2,...,, where n is the number of nuclei in an — 
elementary cell) and by the displacement vector uj,. The origin is placed at a 
nucleus within the crystal and the instantaneous position of any nucleus is then 
specified by rj, =R'+r,-+ ub. 

Separating the motion of the nuclei in plane waves we may write (1) 


h t 
ce? (; Tia) CF (q) [aq ™ 0 — ait he +090) (1.1) 
with G*(—q)=—Ci(q) and «,(q)=0,(—q) (j=1,2,...,32). 


(*) See also: R. J. Glauber, Phys. Rev., 98, 1692 (1955). 
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Here each eigenstate has been labelled by the double index (qj). N is the 
number of unit cells in the crystal and M, denotes the mass of the kth nucleus. 
aq and the Hermitian conjugated one aj, are respectively the annihilation and 
creation operators concerned with the harmonic oscillator. They must fulfill the 
commutation relations below 


Ass —As's=0 | 3=(q/). (1.2) 


+ 
as as. — as As = Oss’ J 


The eigenfrequencies and polarization vectors are determined from a system 
of equations 
wo” Or= > Dry CF, (13) 
ky 
where the dynamical matrix Df* (1) is related to the interatomic forces in the 
erystal and the wave vector q by 


; ; st 
= > ve, (I) 


m (1) being force constants. C* are ortho-normalized as follows: 


kx k 
> Ciq,2 Chaya - Oj7"5 


k,z 


(1.5) 
2. Che fay y = ORK" Oxy. 

The wave vectors are restricted to N discrete values, and if drawn in the 
reciprocal lattice they will uniformly fill a unit cell. Now WN is supposed to be 
very large and therefore a summation over q can be replaced by an integration 
over the reciprocal unit cell. For the same reason it is convenient to introduce 
the density function f(w) defined in the following way: 


f(w)d@=the number of eigenstates in the frequency interval (w, wo +d) 
divided by the total number of eigenstates 3n N. 


As shown by Van Hove [39] and Phillips [40] f(@) is continuous and its deriva- 
tive has a certain number of isolated singularities. 
To each eigenstate corresponds a mean energy 


ho, h 


5 coth (fs), p= 2k, T : (1.6) 


és = 


where 7’ is the temperature of the crystal and k, is Boltzmann’s constant, and 
thus the total vibrational energy per nucleus may be written 


£,=3 | “2 coth (Bw) f(w) da, (1.7) 
0 


where w,, is the maximum vibrational frequency. 


(4) See: M. Born, Reports Progr. Phys., vol. IX, 294 (1943). 
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All discussions of the neutron scattering can be based on the time cor- 
relation function 7M}, (t) = <u (0) uj, (t)>7, where <--->z is the statistical expecta- 
tion value at temperature 7’. Using the commutation relations (1.2) and further 
(1.6) it is a straightforward matter to derive the following expression 


iis <ob Sub yeyie 52 aor ob, es (L8) 
ay : 8 


| 


where q; (t) = coth (f ws) 


cos (wf) ca ; sin (Ws t) 
Ws Ws 


(1.9) 


For crystals with cubic structure the tensor M?., simplifies to a scalar times 
the unit tensor(!) and if further restricting ourselves to those crystals which 
constitute a simple Bravais lattice we may write 


Om 


MW (t)= I | [coth (Bw) cos (wt) +7 sin (wt)] He) dw, (1.10) 
0 


J being the unit tensor.(?) The index k is now unnecessary and therefore 
omitted in the expression above. We note that all properties of the crystal are 
contained in the frequency function f(w) and the maximum frequency @,. 


Chapter 2. Differential scattering cross section 


To account for the scattering of neutrons Fermi’s pseudo-potential [41] is 
introduced to represent the interaction between the neutron and a single nu- 
cleus. Generally we must distinguish between scattering with parallel and anti- 
parallel spins of the neutron and the nucleus. We can account for both possi- 
bilities by introducing a spin-dependent pseudo-potential 

¢ 2 
y= 72% (42a! 8-8) d(r—R), (1.11) 


Here r and s denote the position vector and spin vector of the neutron, the 
latter in units of %, and R and S denote the corresponding quantities for the 
nucleus. m is the neutron mass and (3) 


4 gL tie y 
PE Ss rae 
1 (1.12) 
Me Os a 


(1) We may draw the more general conclusion that the tensor Mx has the same sym- 
metry as the crystal structure. 


(?) The scalar factor will later be denoted by G7 (0. 
(*) See for instance: L. Rosenfeld, Nuclear Forces, Amsterdam 1948, p. 106. 
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where 7 is the spin and a, and a_ are the scattering lengths for the bound 
nucleus referred to parallel and antiparallel spins respectively. 

If the study is restricted to non-magnetic materials (1) the interaction between 
the neutron and the whole crystal is obtained by summing the single nucleus 
potential ([.11) over all nuclei. We must then label quantities depending on the 
kind of nucleus by the index (lk).(2) 


2 
W t Be eae ae ae 
e ge is a 2, 0%, 6(&—Fx), (1.13) 


One = On + 2 ax (Sy 8). 
With Born’s approximation for solving the scattering problem the integration 


over the neutron space coordinates gives the following expression for the dif- 
ferential scattering cross section 


do Es 


> > pn, | (| F| 29) |? 6 (w@ — wy + Wnn,)- (1.14) 


The summation is extended over the final states |n) and the initial states | 9) 
of the crystal, the latter weighted by the statistical factor p,,. The summation 
should also include the spin states of the neutron. Further the following nota- 
tions are used: 


wo=E/h, w)o=E,/h, where E and E, are the energies of the scattered and 
incoming neutron respectively. 

EH, =energy of the final vibrational state. 

E,, =energy of the initial vibrational state. 

nn, = (En — En,)/h. 


The 6-function assures the energy conservation. (| F'| 9) is the transition matrix 
and F is given by 
i > On exp (ix-Ti) 
lk 


and x=k,—k, k, and k are the wave vectors of the incoming and scattered 
neutron. 

The scattering cross section can be simplified by introducing a time-dependent 
pair-distribution function [8, 9]. Since we proceed in a similar way as the 
authors referred to a presentation of the details seems to be redundant. We 
account for both the spin dependence of the scattering and the isotopic disorder 
in the crystal. 

The 6-function in (1.14) is first written as a Fourier integral 


§ (0 0+ Onn) = 5 | exp [t (w—,)t] exp [i (£, —Hy,)t/h] dt 


and further 


(1) The non-magnetic interaction between electrons and neutrons is neglected. 
(2?) Compare the notation of the positions of the nuclei in the foregoing chapter. 
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tH 1tH 
exp [i (E, — En,)t/h] (n| F| m9) = (n| exp (“) Pex exp (= 3 Jin). 


where H denotes the Hamiltonian of the system of nuclei, supposed not to 
contain any spin variable. 


itH —ttH 
exp | F exp i 


may be considered as a Heisenberg observable which can also be written (7) 


)= a 0, exp [ix-ri (t)]. (1.15) 


The differential scattering cross section is now transformed as follows: 


oo 


Es. | exp [i (w — 9) t] dt {5 > Vn, (mo| F* (0) | n) (| F (t)| 9)} 
€dQdw ki, 2x mire & 
Krad r A 
oi, oe | exp [t (w— wy) t] dt (F* (0) F (t)>r. (1.16) 


The completeness of the eigenfunctions has been used to perform the summa- 
tion over |n), including the spin states of the scattered neutrons. The second 
summation is contained in the definition of the statistical average value. As 
shown below this is independent of the spin states of the incoming neutron 
and thus the summation over the initial neutron spin states only gives a factor 
of unity. 

We further note that 


Su 8>r=<(Six°8)">r=0, <(Si.8)" (Sra *8)>2= & Jue Jue + 1) Ow Oxx- 
and if using these relations 
(Ole Orn > = Cerne rer + Oct Jue Ge + 1) Ow Sick (1.17) 
By separating the spin states and space states of the crystal and putting 
the expression for F(t) in (1.15) into (1.16) this can be simplified further. To 


save space we introduce the notation 


Gir = exp [—ix-r} (0)] exp [ix-rh (t)]>r- (1.18) 


We can then write 


lk v 


<F* (0) F (t)>7= 2 Ohrn Orr ner Giese + > oti Fx (Jue + 1) Gite. (1.19) 


(*) The symbol ¢ then has the physical meaning of a time. Cf. L. Van Hove, Phys. Rev., 
95, 249 (1954). 
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For mono-isotopic nuclei we may drop the index J in x, and «;,. When there 
are several isotopes of a certain kind of nucleus we instead define the mean 
values, <a> for instance, over the different isotopes with regard to the abun- 
dance.(?) 

Owing to the assumption of complete isotopic disorder and the equivalence 
of all lattice points(?) it can be shown that 


SB aaetrw Ge NS Cond <a> (5 Oe) +N 5 (<a> aud) OB 


and analogously for the second term in (1.19). By further using (I.12) we get 


ems (0) P(t)yr=N > (Oe) (hie) (> Gv) +N > (<ai> — <ax>*) Gtk, (1.20) 


_ where <a,> and <az> are respectively the mean values of the scattering length 
and its square, averaged over the two spin states and all isotopes. 

The scattering is in this way separated into two terms, the coherent part, 
which includes all interference effects, and the incoherent part, which has the 
character of one-nucleus scattering. 

As shown by Van Hove [9] and Glauber [8] we can finally rewrite G%!,, in 
the following form (°) 


Ghee = exp [—4%- Mx (0)-%—-4b > Mex (0) H+ %* Mie (tx) (1.21) 
Using (1.20) and (1.21) we get from (I.16) the following expression for the 


differential scattering cross section per unit cell, separated into the coherent 
and the incoherent part: 


oo 


> A, exp[-*- Mex (0) xh 52 | exp [i (@-apo)t] exp [> Mix (t)-%] dt, 


2 
d Gincoh _ k 


€dQdw ky t £ 
= (1.22 a) 
@ Coon ko ; —1%- me, (0)-x— 
==): k Lhe? ? 2 
dQdw ky 2, <M) Bx) exp [1% + (Te —¥;)] exp [— 3% Mix (0) +x 
—4- Mew (0)-%] > rent = [ exp [i (mo) t] exp [%- Max (t) >] dt. | 
i J | 
. —oo 
(1.22 b) 


Here A,,= <a; — <a>. 


(4) The mass difference between the isotopes causes a distortion of the plane waves in the 
erystal and therefore we must assume that isotopes with large mass differences are rare. 

(?) The equivalence is fulfilled only if surface effects of the erystal are neglected. 

(3) The proof is based on the assumption that the anharmonic effects on the crystal 


1 : , 
vibrations can be neglected. The tensors Mix were defined in chapter 1. 
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The phonon expansion of the scattering cross section, in which each term 
contains processes where a fixed number of vibrational quanta participate, is 
obtained from (1.22 a, b) if exp (x: Mix (t)-%) is written as a power series. 
The integration over ¢ is then easily performed and the summation over / leads 
to generalized 5-functions, periodic in the reciprocal lattice as 


dp(x)= >, d(x—7). (1.23) 


« is a reciprocal lattice vector defined in such a way that (c-R’) is a multiple 
of 27. dp(x) is connected with the J-summation by 
1 (2a)° 


2, = Sp (xm), (1.24) 


where V means the volume of a unit cell in ordinary space. 
For surveyability we only consider a Bravais lattice and then we get the 
following expression for zero-, one- and two-phonon terms [6]: 


don :, 9 & (22) 
ADomeos?) Bak 


h|x-C,|? coth (6 a,)+¢€ 


‘ exp [—H- T° (0) >] ]3p (90) 8 (0 — 0) + 


2, NM a, 9 On (K+ qs) 6 (@ — Wo + EWs) + 
1 { h|x-C,? h|x-C,.|? coth (8 w,) +e coth (Bos) +e" | 
2 ssvee’ 2NM ow, 2NM ws: 2 aS 


x dp(ux+q,+qs-) d(@—@yteo,t+ é' ws) + |, (1.25) 


é= +1 corresponding to emission and absorption of phonons respectively, and 
summation should be performed over both values. exp (—x-M°(0)-%) is the 
Debye-Waller factor and will also be denoted by the more common notation e?™, 
The incoherent cross section is obtained from (1.25) if unity is substituted for 


(220)° 


Op(...) 


and A for <a*). 


Chapter 3. Total scattering cross section 


Most of the earlier investigations on neutron diffraction were concerned with 
the total scattering cross section. In the first calculation to obtain the tem- 
perature dependence of this cross section Weinstock [5] only considered zero- 
and one-phonon processes. Later Squires [11] derived approximate formulae for 
multi-phonon processes of any order, which he compared with available experi- 
mental results. These show that in many cases a great number of terms in 
the phonon expansion are required to give enough accuracy. In several papers 
Placzek [7] has subsequently suggested that it would be more convenient to 
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expand the total cross section in powers of m/M, m being the neutron mass 
and M the nucleus mass. In this way he obtained simple formulae for the in- 
coherent cross section, applicable to large and small neutron energies. In an- 
other paper Placzek and collaborators [27] also show how to correct for the 
interference effects when dealing with rather high neutron energies. Later Ko- 
thari and Singwi [26] used the m/M expansion for the inference terms as well. 

It is possible to make an exact integration of the differential scattering cross 
sections (1.22 a,b). We need only consider the coherent part; the incoherent 
cross section is then obtained as a special case of the coherent one. 

To abbreviate the notations we write 


R=R'+r,.—r,, 
M=} Mex (0) +4 Mer (0) — Mix: (t). 2) 
Then we get from (1.22 b) 


dean = 3S Ks) <a) sae | exp [i (w — wo) t] dt {e'** exp [—x-M-x] kdQdw}= 
Tk 22k, 


h ; ht 
= 2 2, <a> <aK> oink | dt [exp E . (R Pi. k,)| . 


‘exp | —x:(m—i 2 9) -x dx} 


1 1 r i rf Hein: 
ne Paes | || B || exp [-4(R’-B-R)] dt, 

RR" potnm_-itty 

where ont SY LL tael ye 


and ||B|| means the determinant of B. 
For surveyability we specialize to a Bravais lattice. We then get the follow- 
ing expression for the total scattering cross section [42]: 


: 00 | 
oy hy po : ex Shey “(P= it) *-I,| ———— | (1.27 a) 
Oincoh — ai Pp 2m ° 0 (| — tt J||)*?) 


0 


: i oo i ; ; : oy 
(eon = <a) (2) 2 | ex | - se (FRI ths) rae) ‘(7 R'1,)| : 


| sf Sine 
| a EELS 
(1.27 b) 


(1) We may always consider m symmetric. 
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1 _ 2M ano l 
Here ie mage (M° (0)— M(t). 


At a later stage we will specialize further to cubic lattices and then 


T= dbo { [coth (8 w) (1 — cos (w#)) —% sin (w#)] {ta} da. (1.28) 
0 
The integrands in (1.27) are singular for t=0 when /=0 and some caution 
is therefore necessary. It is possible to avoid the difficulties by putting 


exp [i (@— wt] +> exp [i(w—@,) (t+7e)], &>90 


in the proof and after the t-integration let ¢—0. 


Chapter 4. Discussion 


The derivation of the formulae for the scattering cross sections may appear 
to be very general. Nevertheless there are some limitations which must be 
emphasized. 

Firstly the harmonic approximation of the vibrational motion of the nuclei 
is not justified for high temperatures. It is certainly possible to account for 
the thermal expansion of the crystal, at least in principle, by considering the 
nuclei to vibrate around their equilibrium position at any temperature. The 
interatomic forces are then supposed to be temperature-dependent and, for in- 
stance, the frequency function will vary with the temperature. It should also 
be mentioned that for metals the basis of the theory of the nuclear vibra- 
tions is not all clear. 

The influence of the anharmonic forces in the crystal may be divided into 
two different effects. The first effect manifests itself through an expansion of 
the crystal still keeping the harmonic nature of the motion and the formalism 
above will be unchanged. In consequence of the other effect the plane waves 
are not exactly eigenstates of the crystal and we must suppose an interaction 
between the separate phonons, which thus have finite lifetimes. Such an inter- 
action is also due to other effects as for instance crystal imperfections. Among 
other things this effect prevents the energy-momentum relation for the one- 
phonon scattering (see formula (I.25)) from being exactly fulfilled and causes a 
broadening of the in the ideal case discrete energies of the scattered neutrons. 
Nevertheless it is interesting to know the results in the harmonic approxima- 
tion even when the assumption that the anharmonic forces are negligible is not 
completely justified. 

The other limitation to be mentioned concerns the Bragg scattering. To 
obtain accurate values on the intensity of the Bragg peaks we must, except 
for very small single crystals, account for primary extinction, and this can only 
be done by using the much more complicated dynamical theory of scattering 
[43, 44]. Now it is possible to separate the Bragg scattering from the other 
scattering, both theoretically and experimentally. When dealing with the total 
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scattering cross section polycrystal materials are often used, where extinction 
effects can be neglected. Except for the coherent elastic part the influence of 
the extinction effects on the results should be very slight. 

The one-phonon scattering intensities, both coherent and incoherent, are very 
sensitive to the detailed vibrational properties of the crystal [22-24]. Any rough 
approximation of the eigenvibrations, such as the Debye approximation, would 
therefore fail to exhibit many of the most characteristic features of the one- 
phonon scattering. It is moreover possible to perform calculations on this scat- 
tering by using more realistic assumptions. As shown by the author in another 
paper [34] the coherent two-phonon scattering cross section, too, is rather sen- 
sitive to the details of the vibrational properties. The higher phonon terms, 
however, are much less sensitive to the details, and the interference terms in 
the coherent scattering cross section (J+ 0) make only a rather small contribution 
to the scattering. 

We note that the term corresponding to 1=0 and k=k’ in the coherent cross 
section is equal to the incoherent cross section apart from a numerical factor. 
In this paper we restrict ourselves to this term (incoherent approximation) 
leaving the interference part for a separate investigation. 


Parr II 
DIFFERENTIAL SCATTERING CROSS SECTION 


Chapter 1. The Doppler effect 
Cubic Bravais lattices (1) 


The characteristic features of the scattering are very different for “large and 
small neutron energies”’ (?) and therefore especially the latter of these two extreme 
cases will be discussed before passing on to intermediate energies. 

For “‘small energies” the zero- and one-phonon terms (see formula 1.25) give 
the main part of the scattering, the remaining terms only causing a weak and 
smooth background. In the case of “large energies”, on the other hand, many 
phonon terms contribute appreciably and here the nuclei may be considered free, 
having a prescribed momentum distribution. The deviation of the scattering 
from that of fixed nuclei should thus mainly be due to the Doppler effect.(*) 

The zero- and one-phonon scattering have been thoroughly investigated in 
other papers [20-24] and therefore we need not discuss them further. The other 
extreme case has also been treated before [4, 8] but will be considered in more 
detail. 

(1) The cubic Bravais lattices are the simple cubic, the body-centered cubic and the face- 


centered cubic lattices. 
(2?) This should mean a division into large and small values compared with unity of 
2W= ((x- a which denotes the mean square displacement along % of a nucleus multiplied 


by «2. The neutron energies considered in practice extend from some few eV down to about 


10~* eV. gl 

(*) It should not be assumed that the interference effect on the coherent scattering is 
negligible for large neutron energies. On the contrary, Placzek and collaborators (7, 27] 
have shown that for the total scattering cross section the interference correction will for 


heavy nuclei exceed the Doppler correction. 
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For convenience all frequencies appearing in the following will be measured 
in units of the maximum vibrational frequency ,, and the time in units of 


1/@m. 
hx h 


Thus, for instance, C-omM and P30 (ILAj 


are dimensionless quantities. 
The differential scattering cross section was given in (I.22a) and for a cubic 
Bravais lattice this expression simplifies to 


d” Oincon__4 & a lie [i (w—@ ) #] exp (Q[ (t) meas (0)]) dt (11.2) 
dQdw ~ky2x . P Ki : 
1 
where y w= | [eoth (BC) cos (€t) +7 sin a) iat (IT.3) 
0 


depends on the vibrational properties only by the frequency function f (¢). 
Qy(0)=2W is the ordinary Debye-Waller function and when 


2W>1 


the region near t=0 makes the main contribution to the integral and there- 
fore y(t) may be expanded in powers of ¢, neglecting in the first approximation 
third and higher order terms. 


Ald bt a hah Brees ; Pee, (11.4) 


Ex, = E,/2 is then the mean kinetic energy per nucleus. 
In this approximation, henceforth called the Doppler approximation, the time 
integration in easily performed and gives the result: 


| Porson, & (22)"* (@ — w+ Q) IL5 
\dQdo “% A exp | 2 A? | ie 
4Q Brin 
h Cr! ——ae / 
where A ae (I1.5’) 


For high temperatures Hin approaches the value 3k,7'/2 corresponding to 
the kinetic energy per nucleus for a Maxwellian gas of free nuclei. The scat- 


tering will be approximately the same as for neutrons scattered against such 
a nuclear gas.(1) 


(*) Cf.: B. N. Brockhouse, D. J. Hurst, Phys. Rev. 88, 542 (1952), formula (5). 
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As shown by Brockhouse and Hurst [13] this free gas model gives a fairly 
good interpretation of their experimental results on lead and aluminium at the 
neutron wavelength 0.48 A (H=0.35eV). In (II:5) the model is completed in 
as far as the actual momentum distribution of the nuclei in the crystal is in- 
serted, a fact which is proved in Appendix I.(1) 


General lattices 


In the general case it is not possible to express the scattering cross section 
in the vibrational frequency function. We may, however, expand 7M}, (é) in 
formula (1.22 a) in powers of ¢ and further proceed exactly as above. The cross 
section may be written 


d* aS DS a” ox , 
€dQdwa —~dQdw 


(11.6) 


where each term on the right side gives the scattering from nuclei with a defi- 
nite basic index k. For d*o,/dQdw we get 


d* o;, k (2x)* z ena Oe) IL7 
FOG hy eee a BAR pert 
hx 2 k 
Here Qe = ou ? At =2 Qe (Cy & ey), (11.8) 
k 


where ey means the unit vector along x and ¢«” is the (3x8) tensor (?) 
1 * 
ky ee CE Ob ot (11.9) 


e* has the same symmetry as the crystal and therefore for crystals with 
cubic structure, generally containing several nuclei in each unit cell as, for in- 
stance, NaCl, «* is a diagonal tensor with all its diagonal terms equal to 


1 
K — —__ CGP IT.10 


which after multiplication with ne may be interpreted as the mean kinetic energy 


of a nucleus with the basic index k. 7 
e* can also be expressed by the components of the dynamical matrix [Dr, ] 
as follows (see ref. [45]). 


de-a [aw coth (6 D*) |i dq (11.11) 
cell 


(1) The proof is based on the harmonic approximation of the interatomic forces, so the 
“actual momentum distribution’? means the distribution obtained in that approximation. 
(?) See formula (1.6). 
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and for B<1 and f>1 this expression simplifies to 


: 1 pele | kk a: ble" | 2ykk g mete IL.11’ 
p<]1; eh 5 [ont 3 (22) } ore 45 Qn) J [D dP q iF ( ) 
p>1: peor | wstaq. (11.11”) 


cell 


From (II.11’) we see that for 6<1 the difference between the scattering from 
the separate nuclei is mainly due to the mass difference. “ 

Using (1.4) the q-integrations in (II.11’) can be performed, giving a direct 
relation between the scattering and the interatomic force constants. 


a | Daq=V (0), | 

, i (11.12) 
ee | ae 1) B(V’). 
(220)? jr aq ee aA, 


(i) are then considered as matrixes with the components U%}' (1) and in the 
latter relation these matrixes are multiplied by each other, for which R'+ R'’=0. 


The total scattering cross section 


Though the total scattering cross section is discussed in a,separate part of 
this paper, some results obtained in the Doppler approximation will be given 
in this connection. Having first considered a cubic Bravais lattice we can easily 
generalize the results to other lattices. 

The total cross section given in (1.27 a) is first rewritten in the form 


se eV a [i aot] dt (II.13) 
rm A (5) J re. Saal (oF 
where 0 (t)= —it- Fly ®-y (0)]. . 


To get the Doppler scattering we must expand y (¢) to the second order in t. 
After a simple transformation of the integrand we may write 


jeexp [wA?/(tt+ wA)] dt 


inco =A 3 At ap li 5 
ce ae ee (it+uat  (—it- +8) aan 
; 3 Mw m. | 
with = Md mip ite ot en 
om Ba” aes | (11.15) 
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6 has been introduced to avoid the singularity at t=0, according to the re- 
marks below formula (1.28). 
Using the relation () 


ety 5 izt exp [wA*/(it+ wA)] 1 —pAt .° en 
ak (it+ way! a (nate sinh (24Vux), «>0 pas 
=—() oO 


and Plancherel’s theorem for Fourier transforms we now get 
Gincon = 4A re* yw * Ate * { xt e-** sinh (2AV uz) da, 
0 


and when 2’=2? has been substituted and a partial integration made, this 


formula transforms to (?) 


| 42 A iey 2 1 | 
ei eee ee £ “Al. aly 
oe we It a nt Erf (A =r ey é | | ( ) 


For A>1 and A<1 we may also write 


4A 2 

A<1: oincon=—g- =] (A P+ 4AE+ ++), (I1.17’) 
fe 

42 A il ” 

A>1: A aL ita i (1+55) : (U7 ) 


Comparing (11.5) and (II.7) we can conclude that the above result may be 
generalized to general crystals with cubic symmetry in the following way: () 


Oincon = Cis (11.19) 


where each o; is obtained from (II.17) provided the following substitutions 
are made: 


mm 
AA, pl se eae (11.19) 


For non-cubice crystals wo do not obtain such simple results without making 
further approximations. 


(1) The formula is taken from G. A. Campbell and R. M. Forster, Fourier Integrals, New 
York 1948. : 

(2) Erf(#) means the error function. 

(°) See the remarks on the tensor €* below formula (II.9). 
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Formula (II.17’’) was earlier derived by Placzek [7] and Wick (28] without 
assuming the harmonic approximation, and they also gave the correction due 
to the direct binding between the nuclei. Their results may be written in the 


form 


4n, A m Ean m ye? C | 
aber Leet @ 11.20 
Meier) io f 3 Miw, 32M (hw,) fer 


2 


i 
C is defined from C= 5 p< (gtad V)*>z, 


where V means the potential function for a separate nucleus. In general this 
formula is valid for w)>1 except for very light nuclei and high temperatures, 
where the complete expression for the Doppler scattering (II.17) should be used. 


Chapter 2. The phonon expansion 


Unsymmetric formulation 


If the neutron energies and the scattering directions are such that 2W is 
of the order of magnitude of unity neither the higher order nor the low order 
phonon scattering may be neglected. The method to be used must be developed 
in such a way that significant features of the scattering which depend on the 
detailed properties of the crystal vibrations are not smoothed out in the ap- 
proximations used. It is then convenient to take the phonon expansion men- 
tioned before as a starting point for our treatment. The expression for the scattering 
cross section is first written in a form suitable for direct numerical computa- 
tion and from which we can also see how to make simplifying approximations. 
Two different formulations, called the unsymmetric and symmetric ones, are 
given. Whether the first or the second is most convenient for use depends on 
the temperature and the scattered neutron energies considered. At first we re- 
strict ourselves to cubic Bravais lattices, making generalization to general lat- 
tices in a later chapter. 

The phonon expansion is obtained from (II.2) by writing 


2 
exp (Q[y (4) — y (0)]) =exp [— Qy (0)] j1+ey w+ cy (t) + -.| 


and determining the Fourier transform of each power of y(t). Thus the trans- 
form of y (t)/y (0) 


9 (0 09) = == J exp [t (@ — wp) ett (11.21) 


may be written 
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g (w) ar ? | @ | S i. 
oy (0) 2 | (1.22) 
| = 3 60) Si! | 
f(w) is then defined for negative frequencies by 
f(—w)=f (a). (11.23) 


We notice that g(w) has the same kind of singularities as f(w) and in gen- 
eral shows pronounced maxima and minima(!), and that g(w) is normalized 
as follows: 

[ 9(w)dw=1. (17.24) 


—-ce 


The Fourier transform of any power of y(t) can be expressed in g(w). We 
may write 


— 
~ 
— 


1 ? : i 
Gn (@ — Wo) = 5 | ex [t(@— @p) t] Ral dt = 


= | g@-eng oo, soe GALn ad — Wy) OO, «20 nt 


The limit of integration is only formally extended to infinity, for the integrand 
differs from zero only within a finite interval. We notice that 


G,(w)=0 for |@|>n and | G,(w)dw=1. (11.25) 


The scattering cross section may now be written in the following form: 


| 72 ro) Ww)" 
| a Oincoh k —2W (2 ) 
Oe G = A IT.26 
dQdo See aah Oe rei 


G,,(w), which are independent of the scattering direction, are determined from 


(4) See the frequency functions given by Walker [31] and Squires [23] for Al and by 
Jakobsen [32] for Cu. 
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Go (w) =06(a), | 

G, (@) g (a), 

G, (a) = | g (@— @ ) Gy (w ) dw’, (11.27) 
Gnii(o)= | g(w—o') G, (w') do’ 


This formulation is unsymmetric in as far as G,(—@)+G,(m). For high 
temperatures the asymmetry is small but at the absolute zero G,(w)=0 for 
w>0, owing to the fact that phonon absorption is then impossible. 

The first term in (II.26) corresponds to pure elastic scattering and the in- 
tensity depends on the crystal vibrations only by the Debye-Waller factor. 
The second term is very simply related to the vibrational frequency function 
and, as especially emphasized by Placzek and Van Hove [22], this fact can be 
used to determine the frequency function experimentally by measuring the 
energy spectrum of the neutrons scattered in a certain direction. Such experi- 
ments have recently been performed with fairly good results [15, 17-18]. The 
higher phonon terms are then depressed by using cold neutrons. 

The multi-phonon terms are, however, much less sensitive to detailed assump- 
tions regarding the crystal vibrations and it would therefore be possible to per- 
form fairly accurate calculations of these terms without detailed information 
about the crystal properties. The Debye model would certainly be accurate 
enough for most purposes but it does not greatly reduce the numerical work. 
In the next chapter we will therefore use another method, which has been 
applied in mathematical statistics to calculate Faltung-integrals as appearing in 
(11.27). 

For 2W<l1 the phonon expansion converges rapidly but for larger values 
of 2W the convergence is slow and a method of summing up the main part of 
the phonon terms will be given. 


Symmetric formulation 


When (II.26) is used for direct numerical computation it is often more con- 
venient to transform this expression as follows to get a more symmetric form. — 
We write 


= pha f() 
ios ¢ 2@y (0) sinh (Ba) 
g (@) 


and transform (II.26) to 


Gone aoe (2 W)" 
SO A aaa 10g) 
| dQdo ae Parr 


G, (w — a). (1.28) 


G,(@) is determined analogous to (11.27). 
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f(@) 


|w|<1, 


=(), |w|>1, (11.29) 


Gnai(@)= | 9 (@— a4) Gy (o,) doy. 


—-o 


In that case G,(—w)=G,(w). The whole asymmetry appears in the factor 
exp [—B (o— a) ]-() 

Especially for high temperatures, when the asymmetry is small, and further 
if we wish to determine the intensity of neutrons that have gained energy in 
the scattering it seems more convenient to base approximations on (II.28) rather 
than on (II.26). 


The Debye approximation 


As mentioned before the Debye model for the crystal vibrations is certainly 
a good approximation for the higher phonon terms, and for high temperatures 


further approximations can be made, after which G, can be easily calculated. 
g(m) is expanded in powers of /. 


3 
2B y (0) 


and for 6<1 we may thus approximate to 


(w) = 


QI 


ae 
(1-F ots ~), |w|<1 (11.30) 


(11.31) 


y (0) has then in an analogous way been approximated 


y (0) = ; = = , ©= Debye temperature, 


If we now write G@,(w) in the following form, which corresponds to the 
approximation in (II.31), 


= Lee sin t\” 
G,, (w) ==— J ere (=) dt, (1.32) 


the integration is easily performed by calculus of residue. We get the fol- 
lowing result: 
(4) Note that G,, (@) is not normalized to unity as G,(@). 
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((n—|o|)/2] (n—2y— | wo ie 
bch < 
Pal yLiwelieo Sa 


The summation is extended over all integers from zero to the largest one below 
(n—|o0))/2. : 

We notice that in this approximation G,(w) are independent of the tem- 
perature. The figure above (Fig. 1), where some of the G, are drawn, shows 
the functions for n>2 to be very similar to Gaussian ones. This fact will be 
used in the next chapter to calculate the high-phonon scattering cross sections 


for general lattices. 


Chapter 3. Expansion of G,(@) in Hermite polynomials 


Unsymmetric formulation 


It is a well-known fact that G,(w) always tends to a Gaussian function as 
m increases.(1) The same is of course also valid for G,(qw), as illustrated in 
iigs 1: 

To account for the deviation from a Gaussian function G,(w) is expanded 
into Hermite orthogonal functions as shown in the following. The method is 
described in textbooks on mathematical statistics and therefore a brief sum- 
mary is sufficient in the present context.(?) Owing to the fact that G,(w) is 
different from zero only within a finite interval we have no trouble with the 
convergence of the expansion. The essential thing is, however, that even in 
very accurate calculations we need not consider more than the first few terms 


in the expansion and for more approximate calculations all but the first term 
may be neglected.(3) 


(*) This theorem is known as “The central limit theorem’. See for instance: H. Cramér, 
Mathematical methods of statistics, Uppsala 1945, p. 213. 


(2) See for instance: H. Cramér, loc. cit. p. 185, 221. See also: C. G. Essén, Acta Math. 
77, 1 (1944). : 


(*?) We only treat G, for n>2 in this approximate way. 
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To get a convenient starting point for our treatment G,,(w) is written in the 
following form (see the foregoing chapter): 


= ; ic (t) dt 11.34 
ees e [2 ale : (11.34) 


We further notice that ly ()|<|y(0)| for +0, (11.35) 
except for an ideal Einstein crystal, where y (t) is strictly periodic in t, and that 
y(t)=o(€") when t>co. (I1.35’) 


The latter property is due to the fact that the vibrational frequency function 
is continuous and its derivative is integrable.(1) For large values of » a small 
region around t=0 makes the main contribution to the integral except possibly 
for large values on w, where, however, the function is very small. y (t) is there- 
fore expanded in powers of t. We may write as follows: 


O iets a (ity, (IL.36) 


where the coefficients «, are determined from the vibrational frequency func- 
tion by the following relations 


| fs 
—— 


=I, fare (0), 
1 . 1 
pei _ 4 Brin Plz bes - 
1-5 [ecomipetcrac= 5 id ma Heaes| arse 
0 
1 
end 
OL coth (BE) f(C) dé a scaag te 
The corresponding expansion for In [y (t)/y(0)] is written 
YQ) | _ S Heyy 
In 2 2 | 2 51 (oe) (11.37) 


and x, are then determined from the «, by 


My = Oy — Ot, 

11.37’ 
Hg = yg — 3. % + 2a, ( ) 
Hq = 04 — 30R— 40, %3 + 12 af a — 6 act 


ww Oe Oe ae eS re Oe ew [et a  @ e”8 


(1) Van Hove [9] has shown on rather general assumptions that y (t) = 0 (é~ 4). = an ideal 
Debye crystal the frequency function is discontinuous and in that case y (t)=0(t7}). 
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A : . 2 ree 
For convenience we introduce a special notation Ao for x, and this is thus 


defined by 


(11.38) 


A [zee 1 |’ 
8 [3h 0) fy OPI | 


G,(w) can now be treated as follows: 


n 
22 ee 
1 point 70) = *( it J] 
eat eg ao =~ 6) @xps| Ri Do; =) |dt= 
on A, Va a ae aN eS ele haiet Laesich GF 


o) 
1 (n) 


~ oak, Va et : 


1 
alk. V2nn 


0 


i x. 
aa f 1 ——+i1—*~t)|dt= 
| ie Sloe aves as 


-o 


exp (— 2%/2n AG) 5 (—1)’ co H, (a,/Ao Vn). 


v=0 


% denotes w+n/y(0). H,(x) are the Hermite polynomials.(!) The first few in- 
cluded in the series above are given explicitly below. 


The cS” coefficients are 


(11.39) 


H, (x)= 2° —15a*+ 452-15, 


determined from 


(*) The Hermite polynomials are defined in agreement with the definition used by Cramér 


(loc. cit.) and also found in 


d”™ 
Tabl 2)—* — (e-T*/2 
‘ables of (27) ant (e ) 
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2/2 d” 
e— 7? H(z) =(— 1)" dat fe“). 


may be found e.g. in Cramér’s book, p. 557. 
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n 
=1, cP=cS=0, 


Ces, ee ee 

3B!AR nt nt’ 
(nm) ee 1 vig: 

4!Agn n’ 

(11.40) 

(n) __ 5 i 
abuts LAR nh? 
tits Pe | “a 


2 “61 AS m2? BIAS? 


PKS) cor hler ste (el) Wet tee. Ade de he 


The symbols y, and y, are introduced to get a shorter notation. They depend 
on the temperature and the frequency function and can be tabulated provided 
the latter is known. 

If the terms in the expansion above are arranged in powers of n=? we get 
the following expression for G,(w@) including only the first three terms: 


| prepa. lie 3 (— 22 /2n AB) ){1— Fa (2%,/Ao Vn) + | 
| a (11.41) 


+5 (aH (tu/o Hn) + 7H He 4/Mo Vm). 


As already mentioned x, is given by z,=a+n/y (0) 

The first term in (II.41) will in all applications give the main part of the 
scattermg and in many cases we may neglect all other terms. Owing to this 
fact it might in almost all cases be sufficient to use the Debye approximation 
when computing y, and y, as long as one keeps in mind that uncertainties in 
y(0) and A, will have quite a great effect on the Gaussian function. In that 
case y, and y, depend only on f# and may be tabulated once and for all. 
Tables of them are found in Appendix II. 

It may, however, be essential to have more exact values on y(0) and A, 
and in calculating these more realistic assumptions on the frequency function 
should then be made. For f<1 and for f>1 we may write as follows, if (w”> 
denotes 


1 
<o") = J wo" f(w)da: (11.42) 


(o> B_ BF 
y(0)= B ts 1p t” ¢ 
aia", B>1, 
-2 2 
[y (0) Aq?= Se £ (498) <a?) —2) +B (A¢0%) — (o> Coty) +o, P<, 
Ss aT p>l. 


ae | 


(11.43) 
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From this expression we conclude that especially for high temperatures it is 
possible to use the Debye approximation also in determining y(0) and A, pro- 
vided that the Debye temperature is conveniently defined. To get good values 
of both y(0) and A, we define the Debye temperature as follows : 


hop nV3 


OT mie 


(11.44) 


In chapter 1, however, we saw that for high neutron energies y (0) disappears 
from the expression for the scattering cross section. Our results will thus be 
quite insensitive to the values adopted for y(0), provided that the same value 
is used in the Debye-Waller factor. In that energy region it would therefore 
be more accurate to define © as follows, giving a more precise value for Ag. 


2 
@= nV <w?)! (11.45) 
kp V3 
Tables of y(0) and A, obtained in the Debye approximation are given in Ap- 
pendix II. 

In a later chapter the approximation (II.41) will be compared with exact 
results obtained from (II.27) by direct numerical integration, assuming a certain 
form of the vibrational frequency function. 

Let us, however, first return to the first term, the Gaussian term, in our 
expansion above and show how to sum over » to get an expression which is 
more convenient for the calculation of the main part of the two- and higher 
phonon scattering. 


The Gaussian term 


We would have obtained the Gaussian term without corrections if we had 
in (II.27) simply changed g(w) to 


(0)= api ow |- (0+ |" /2.03 
I A, (22) xp o ey i. (IT.46) 


The actual function g(w) is then approximated by a Gaussian function with the 
same mean value and the same spread as the original function. 


oo 


<w> = | og(w)do= “70” <(w — <)> = i (@ — <>)" g (w) dw = Ao, 


-o 


J9(o)do= fq (o)do=1. 


338 


ARKIV FOR FysIK. Bd 14 nr 21 


If the Gaussian terms are inserted into (II.26) the sum over n may be 
transformed as follows: - 


2 (2W)” 1 erent ON ce! on bit S z” exp (—2"/2n) 
naa 1! Ag(2an)* AG neet) (2an)* ~ 


z, « and f,, respectively, are defined by 


2 1 1 


2-000) exp (— 3/31), t= 


These quantities depend on the temperature since y(0) and A, depend on 
temperature. The above sum defines a function F(z,x) which is evidently in- 
dependent of the structure of the crystal and may thus be tabulated once and 
for all. In appendix II this function is given for a series of values of z and z. 


(11.48) 


The scattering cross section may now be written in the following form: 


a” Cin h k Qy(0) 1 
con __ ee a ~S — cee aii 
SOds ote Ja ®))+Qy(0)g(@ Mo) aie Fie; 2) + Foor 


(11.49) 


Forr Genotes the correction due to all but the first term in the expansion 
(11.41). Where only an approximate calculation is required we may neglect this 
correction term. 

In such approximate calculations we may also for temperatures higher than 
the Debye temperature make further simplifications as follows: 


67 1 _O¥3 
-—. 


y(0)= 6? Ay= Va By (11.50) 


In determining z we may put /,=0. 

It will finally be proved that for large Qy(0) the Gaussian approximation 
correctly gives the Doppler scattering earlier discussed. To show this the Gaus- 
sian terms are written in the form 


1 2 2 1 | : 2 42 
ee a = — t— t") dt 
A @unt exp (—27;,/2n Ao) 5, ) °xP (ia,t—in Ajt*) 


and subsequently inserted into the formula (II.26). Owing to the fact that the 
one-phonon scattering decreases rapidly for increasing values of the Debye-Waller 
function we may also treat this term in the Gaussian approximation without 
appreciable error for the whole scattering cross section. We then write as follows. 
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oo 


exp [2 (w— @p) ya | es “i J27700) exp ley —ase) |" f = 


(w-@ 
bone 


~ (2W) 
ee n! 


no! 
on 


$0) oe 
exp [7 (w—,) t] exp jor (0) exp (7, -548¢) | dt 


J 
J 


If we now proceed in the same way as in chapter 1, expanding 


Ae 2 42 Lie PN hoe )e 
exp es ms e)- 1+450) 5 (43+ on ie 


we get the folloving result after the integration over f¢: 


W exp [—(w—@,+ Q)?/2A?). (11.51) 


(@ — Wo) = 


1 2 
= A (2a)* 


A is then the same quantity as defined in (II.5’) and the formula above cor- 
responds exactly to that in (II.5). We see how the factor e?” cancels the 
Debye-Waller factor and how the quantity y(0), being essential for low ener- 
gies, disappears in that case. 


Symmetric formulation 


As mentioned before it is, where possible, more convenient to base the ap- 
proximate treatment on the symmetric formula (II.28) and we will therefore 
discuss this case as well. We may proceed in exactly the same way as above 
and need thus only give a very brief presentation emphasizing the modifica- 
tions to be made. 

g(w) is first written in the form 


oP 27 0) i we (200) 
g (w) a )2a J ° 5 (0) dt, (11.52) 
f(¢ 
res (Cc) , 
where p(t)= l Esinh (BO) cos (Ct) d¢. (I1.62’) @ 
Analogous to (II.34) we have 
ata ne | ta pa 
(w) E (0)| 2x » e 5 (0) dt. (11.53) 
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where all odd coefficients %, %,... now vanish and the even ones are deter- 
mined from 
1 

J Al ¢f(c) 

x, = 1; & = - | : the 

> (0) J sinh (86) ¢ 

i “ (11.54) 
1 3 
= $ He) EC nna wigtst ore a. Fae 2%, 


7 (0) J sinh (8¢) 


The coefficients corresponding to x, are determined from (II.37’) with %, in- 


serted instead of the «,. The quantity A >, corresponding to A,, is in that case 
defined as follows: 


A} = & = (2 dé. (11.55) 
0 


The expansion of G,(@) in the Hermite polynomials may be written as 


Gn (@) = 2 eo w*/2n A?) 5 (—1) & H, (w/Ay 11.56 
() Be SSR Ur eed (/Mo Vn). (11-56) 


This formula is simpler than the corresponding unsymmetric one in two respects. 
Firstly only &” with even y» are nonvanishing and secondly X,, introduced 
earlier, is in this case independent of n. The even é&” are given by 


A =1, &=0 
Geel Alt ga 


and the formula corresponding to (II.41) takes the form 


al raiders 2/2 As) | PH, A, IF 11.58 
G,, (w) Ee er eee /2n 1+ (w/Ag Vn ( ) 


Inserting this expression for G,(w) into (II.28) the scattering cross section 
may be written in the following form: 


een 4 E exp [- 0710) Aloe) [>(o-0) +2 70)7(0—o9) + 
g r | (11.59) 
| eae. mth 1(%, 4) « | 
z, and 2, are defined by - 
=Q7 (0), ty = (W — Wo) /Ag (11.60) 


and F(z,x) is the function introduced in (II.48). F,(z, 2) is defined in an 


analogous way by - 
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n 


ee 2 
le lea Aue ee F(a /Va) (11.61) 


and refers to the second term in (II.58). 

It must be noted that the formula (II.59) cannot in general be used for low 
temperatures. For in that case the unsymmetric factor exp (—B(@—@p)) has 
such an effect that exp (—£(w—p))[...] has its maximum value for an @ out- 
side the region of validity of (II.58). The formula should therefore mainly be 
used for rather high temperatures. It is immediately seen that (II.59) and 
(11.49) are then approximately equal. The factor exp (— 2), for instance, tends 
to exp (—f(w@—«,)) as the temperature increases. It may also be shown that 
the Gaussian approximation does not exactly give the Doppler scattering for 
2W>1 but only an approximate expression valid for temperatures that are 
not too low. 


General lattices 


The same method as above may be used to derive approximate expressions 
for the scattering cross section for general lattices. We must then, however, 
start from the general formula (1.22 a) and make the expansion 


io) 


exp (n° Mix-*x)= > (4: Me +)" 


n=0 


Subsequently (e,,- 7x (t)- e,)/(e,: Mx (0)-e,) is expanded in powers of ¢ in the 
same way as earlier for y(t)/y(0). e, denotes the unit vector along x. The 
treatment is the same as for simple cubic crystals and therefore we will give 
the results only, restricting ourselves to the unsymmetric formulation, which 
has a larger range of validity. We can, however, easily derive an expression 
aero: to (II.58). Moreover the modification mainly consists of exchanging 
coth (... for 1/sinh (...) in the symmetrical formulation and putting all odd 
creas a below ee to zero. 

In the unsymmetric formulation we get the following expression for the «*, 
nue es a to the earlier coefficients «, and here refer to a certain basic 
index k, 


1 

o=1, af = : 
: *  y4 (0) 
ie 1 af Kk 2 2 k 
a2 Vi (0) V 2, @; coth (Bas Je Cc | Vk (0) (€,,°€ *,),(5) 
ors lel oo? a (11.62) 
= Vx (O) yw ee 5 |e: A P, 
og = : 205 oth (8 w,) |e, CF |? 

Vx (0) Nits S 


(1) €* was defined in (II.9). 
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y«(0) is that part of the Debye-Waller function 2W, which is dependent on 
_temperature and is defined as follows: 


1 
Ve (=> > w;* coth (6 w,) |e, -C*|?. (11.63) 


With» these notations and the ones corresponding to x,, Ag, y, and y, we 
may write the scattering cross section in the following form: 


a” Sincoh _ S F ofron 
dO dau + dO da’ | 
? ofnon k eo W,)" | 
| nco = pa (k) 

(EOS a lel ror ac 


Gs? and G{® are defined by (1) 


coth (Ba,) +¢ (11.65) 


1 1 
GY = BS) 
* (0) > 20s 


Je, CS]? 4 ne 
and G{? are for n>2 approximately given by 


1 
Be =A Gan)t exp (—a%,/2n AG, x) j1— 2s Hy, (an/No,4V'n) + 
0K (2 


i 
+= (a4 Ha (tn/Mo,x Vn) + 4 V3 x He (an/No, x vn) . (1.66) 


Analogous to the earlier notation x, is defined by 


t= 0 Oy + (11.67) 


Chapter 4. Numerical comparisons 


In order to gauge the accuracy of the approximations discussed in the last 
chapter some numerical calculations have been made assuming a certain form 
of the vibrational frequency function. In the first case (11.31) was supposed 
exactly valid, where formula (II.33) could be used to determine G, (w) precisely. 
These were also calculated from the approximate formula (II.41), referred to 
as Approx. A, and in the Gaussian approximation, where the scattering con- 
stants y(0), Ay, y, and y, were determined in agreement with the adopted 
frequency distribution. This case corresponds fairly well to a Debye-like crystal 
in the case of high temperatures. 


co 
(4) J a dw=1 only if e, is held fix. 


-o 
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In the next two cases we consider a realistic frequency function far from a 
Debye distribution, supposing in one case room temperature and in the other 
a rather low temperature. The vibrational frequency distribution adopted is 
that which Walker [31] obtained for Al from measurements on diffuse X-ray 
diffraction.(!) The calculations were performed exactly by using formula (II.27) (?) 
and further in Approx. A and in the Gaussian approximation. The scattering 
constants were determined in accordance with the adopted frequency function, 
and thus the discrepancies between the different calculations are due only to 
the truncation of the expansion in Hermite functions. 


Case I 
According to (II.31) g(m) is supposed to be equal to 
wis Pi goto 
NO fesoa seg Seeio. |o|<1, (11.68) 
=V,; ee |w|>1. 


g(w) is then correctly normalized. The corresponding expression for y (¢) is 


_ 3 sinh (6 +t) 
=, ati tee (11.69) 


and from this we easily obtain the coefficients «,. The calculation was made 
for 6=0.5, and the scattering constants then have the following values: 


1 
5) 7 O84, Ao 0.563, y= — 0.058, y,= — 0.043. 


The results of the different computations are presented in Figs. 2-3. 

(*) Walker’s frequency distribution is adopted only to have a realistic frequency spectrum 
for the numerical calculations. 

(?) The exact calculations were performed on the electronic computer “Besk” in Stock- 


holm. The author is very grateful to Miss G. Hambraeus and Mr. T. Berggren for valuable 
help in making a program for these calculations. 


344 


ARKIV FOR FysIK. Bd 14 nr 21 


G;(@) 


Fig. 3. 
Case II. 
6=0.748 corresponds to the supposed temperature 7’'=300°K and to the 
maximum vibrational frequency w,=5.91x10" sec’. The scattering constants 
have the values 


1 
——_ = 0,163, A,=0.457, y,=—0.062, y,= —0.087. 


y (9) 
f(w)/m*? and g(w) are shown in the first two figures (Figs. 4-5) and Figs. 6-8 
present the numerical results for G,, G; and G,. G, is also calculated assuming 
the ideal Debye frequency distribution and making the computations according 
to the exact formula (II.27). The Debye frequency wp is defined so as to give 
the correct value for Qy(0) as follows: 


y (0)=yp(0)/wp, (11.70) 


where yp (0) is taken from the table in Appendix II for Bp =hwp/2k, T = 0.748 wp. 
From this definition we get wp=0.83 (fp=0.62). The values for the scattering 
constants above, obtained in the Debye approximation, are (+) 


A,=0.455, y,=—0.071, y= — 0.038. 
We see that these agree surprisingly well with the exact ones. 


Case III 
In this case we assume 6=3.00 (7’=74.8°K) and the corresponding values 
for the scattering constants are 


(1) As A, depends on the unit in which the frequencies are measured we must relate A, 


and A, as follows: 
(@ — y)/@,, Ap = (@ — @y)/@p Ap, 


where Ap is obtained from the tables in Appendix II for 8, = 0.62. y,; and y, on the contrary 
do not depend on the units used. 


345 


A, SJOLANDER, Multi-phonon processes in neutron scattering by crystals 


f(@) 
m2 
0 nt 0.5 1.0 w 
Fig. 4. 
—1.0 —0.5 0 0.5 1.0 w 
G,(@) 
Fig. 5. 


a = 0.432, A,=0.320, y,=—0.140, y,=0.034. 


If defining the Debye frequency as in the preceding case we get in the Debye 
approximation (wp = 0.836, Bp = 2.50). 


A,=0.314, y,=—0.171, Y= 0.030. 
g(w) is drawn in Fig. 9 and Figs. 10-12 show G@ 
346 
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Gaussian approx. cae 
a 


Debye approx. 


Some conclusions may be drawn from these numerical comparisons. Thus the 
approximation given in (II.41) seems to give very accurate values except for 
the two-phonon term. The more simple Gaussian approximation also fits fairly 
well and should be accurate enough in many applications. To obtain the two- 


25: 4 347 


A, SJOLANDER, Multi-phonon processes in neutron scattering by crystals 


\\s~ Gaussian approx. 


—3.0 —2.0 —1.0 0 1.0 2.0 w 
G,(@) 
Fig. 8. 
—1.0 —0.5 0 0.5 1.0 w 
G,(@) 
Fig. 9. 
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approx. Ax/-\\ 


—3.0. —2.0 —1.0 0 1.0 w 


phonon term more accurately we must use a realistic frequency function and 
the exact formula (II.27), which moreover in that case is nearly as simple as 
the other expressions given. The use of a Debye frequency distribution does 
not seem to give much better results than our other approximations. 
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Part III 
TOTAL SCATTERING CROSS SECTION 


Chapter 1. Expansion of the total scattering cross section in powers of m/M 
General formulation 


As already mentioned the total cross section can be brought into a con- 
venient form for direct numerical calculations if it is expanded in powers of 
m/M. Except for very light nuclei and high temperatures this series conver- 
ges rapidly and in general only two or three terms need be considered. 
Placzek restricts himself to a Debye crystal (1) and in this chapter some of his 
results will therefore be generalized to real crystals to get an idea of how the 
details of the vibrational properties manifest themselves in the total cross sec- 
tion. In order not to complicate the mathematical treatment we specialize to 
cubic Bravais lattices, where all effects of the lattice vibrations are determined 
by the frequency function. We can, however, easily modify the results to make 
them valid for general crystals. . 

An expression for the incoherent total cross section was given in (II.13) and 
after introducing a function ® (t), defined by 


y (t)—y (0) 


e=(-in +o), o=— Pr’, (III.1) 


the cross section can be rewritten as follows: 


(1) For large neutron energies, however, Placzek gives results valid under more general 
conditions [7]. 
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} co 
onan = A (=| [ion exp (—iwpt) F(t) dt, (IIT.2) 
0. 


a : ® 
where F (t)=(1+ 0)? exp ( ont 5) . (I11.3) 


As ® is proportional to m/M, the m/M expansion is obtained by simply 
expanding F(t) in powers of © and making the integration term by term. If 
{®|<1 everywhere, corresponding to the condition 


m 


67 
au? (0)S1 (or = <1, ©=the Debye temperature), (IIT.4) 


MO 


the series converges for arbitrary neutron energies h wp. 
We thus write 


F (t)=1+(t@yt—3)O+---. 
The final integration is easily performed if it is noted that 


gnt2 at grtt 


tO) 


* a —(n+ 3) ,~izt = 
au [i oD wamestie he Tainaro mend)? (TIL.5) 


=0, x< 0, 


and the total cross section may be put in the following form 


Oincon = 470 A 2, (77) Ons (III.6) 


where the first two o, are given by (1) 
ern = if 


f (@) 


@ 


4 
[eoth (Bw) —1] (1 + 2) ( - 5. do—20,y(0).{ GL 
Wo 2 


Wo 


Large and small neutron energies 


A more convenient form for numerical calculations is obtained by further 
expanding (III.7) in powers of 1/w,(w)>1). The treatment is straightforward 
and the result for o, can be written in the form 


() Cf.: D. A. Kleinman, Thesis, Brown University, 1951. 
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| 


0, ='>, be Jo (w) [coth (Bw) —1]da, (IIT.8) 
n=0 Wo : 
fficient rel 422) (*) (1) 
where the coefficients Cn sla ah 
co=1, =$ &=0, =—e c=» = — se 


A corresponding formula is obtained for wy<1 by expanding in powers of 
w}. We then get 


| o= Ste; wg +03 Wy t+c3’ wet (III.9) 
ee: 

1 
Pi 
Co eae | w' f (w) [eoth (Bw) —1]da, 

= 0 

1 
cy’ = Jor f (w) [coth (Bw) —1] do, 

1 (III.9’) 
og = — 2 for (w) coth (6B w) da, 


0 
1 


ce =% | w * f(w) [coth (Bw) —1]da, 


0 


The linear term in @, comes from the elastic part of the scattering. 


The Debye approximation 


If Debye’s approximation for the vibrational frequency distribution is used, 
a simplified expression for o, is obtained also by making an expansion in powers 
of 6 (B=h/2kgT=O0/2T). We then get a formula which may be used for 


arbitrary neutron energies. We write 


coth (Gola ee ae 


pore (111.10) 
ebay Ne iio 
LMCI 75 : 


The q@-integration is now easily performed and gives 


(1) (es) are the binomial coefficients. 
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| : 
\o4= a LO are HR o) +B Qo ( (og +008. (IIT.11) 


Q, (@») are determined from 


3 3 
Qo (Wo) = oy| (1+ =) fo(1+ 0) + (1-2) fa(1—)~6) | 
need 
1\2 Tae 
Q (00) = (2 2) fy (1+ 9) — (1 -=) f,(1— on)| (11.11’) 


3 $ 
Qo (Wo) =3 @ o[(1+ fa(1-+ 0) + (1-2) fa(1—o)—2] 5 


where fy (x) =22-1 
| (111.11) 
1 


) 
in x)= a — 8 3? + 12 i 
For w)<1 the terms (1—1/a,) “se (1—a,) should be omitted. 


Higher order terms 


To derive expressions for the higher order terms in the m/M-expansion it 
would be easier to start from formula (II.26) (1) and expand the differential 
scattering cross section in powers of m/M before making the integration over 
the frequencies of the scattered neutrons.(?) The integration over the scattering 
angles leads to a set of functions 


®,, (z) = (14+ a)?"*?—(1—a)?"*” (IIT.12) 


and after the integration over w we get the following expression: 


| et 4 (11.13) 
nana ar 7(0)] Be | 
G=1 
on 1 oo| | &,()/2) 4 on) do 16| 
0 
5 te coe fe (V2) {2 G, (w — @) — G («—an)} der + 32| 
O2 199 2 Ws 2%, 0 1 0 > 
0 
Baers 8 ig V“) ae — Wo) — 3G, ( )+4G, (@— a )} dw —18| - 
gar 3? 3 Wo. ¢ 3 (@ — Wo) Ww — Wo 0) 
0 


(IIT.13’) 


(1) To get an idea of the convergence of the series it was, however, convenient to base 
the derivation on formula (II.13) as done above. 
(?) Cf.: G. Plaezek, Phys. Rev., 105, 1240 (1957). 
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For large neutron energies (w)>1) a modification of the m/M expansion can 
be made, as shown by Wick [28], so that the first term in the series gives the 
high energy limit with the exact mass correction, 4 A/(1+m/M )*. We can 
easily derive this expansion from (II.13) by expanding the integrand in powers 
of (1+m/M)t/w, and then integrate term by term. This series represents an 
asymptotic expansion in 1/m, and the first three terms are found in (11.20). 


Chapter 2. Phonon expansion 


General formulation 


As found in the last chapter the use of the m/M-expansion is restricted to 
cases where the factor (m/M)y (0) is not too large. In order to derive approxi- 
mate formulae that do not have such restrictions on their applicability we will 
discuss the total cross section also from another point of view. We then spe- 
cialize to very slow neutrons, a case of great interest, and base the treatment 
on the phonon expansion and the approximation (11.58) for the separate phonon 
terms.(1) In some respects our results will thus complete those obtained from 
the m/M expansion. As (II.58) gives a good approximation for two and higher 
phonon terms only, the total scattering cross sections from zero- and one-phonon 
scattering should be calculated separately, using more accurate approximations. 

To get the total cross section from (II.28) we must first make an integra- 
tion over the scattering angles, which, when expanding in powers of ap, leads 
to the following expression (2): 


: : 2n® : 
fe W)"e Vdo-42|%(0)| exp = 3770) 0| loro, tort 


2 
oa ayy (0) ots E (0) om] 4 | -  (IIT.14) 


G,(@—p), taken from (II.58), is also expanded in powers of w, accounting 
only for the zero-order and the linear term in wy, and the final integration 
over w then leads to Weber—Hermite functions as follows: (3) 


2 = 
fom exp ( -awo—-— =| do =(AgVn)"*#T (+8) e"*"* D_ nay) (Vm) + (IIT.15) 


2n Ao 


a and x are related to the vibrational properties of the crystal by 


a=F7(0)+B, v= (37 (0)+8) Ao: (III.16) 


(7) Note that for small values of w, and not too low temperatures (f ®,) <1) approximation 
(11.58) is more accurate than (II.41) and in addition gives expressions more simple to handle. 
(*) We have 2 W=(m/M)y (0) (0 +@)—2Vmw, cos #) where ® is the scattering angle. 


(3) See for instance: Erdelyi, Magnus, Oberhettinger, and Tricomi, Higher Transcendental 
Functions, vol. 2, pp. 116-130, New York 1953. 
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The contribution to the total cross section by the nth order term in the 
phonon expansion is written in the following form (cf. formulae (III.6) and 
(III.9): 


4nA 
ofncon=—— 4 [of + 9 of +--+]. (111.17) 
0 


The term o{” is further divided into two parts o{"} and o{") where the first 
part refers to the Gaussian approximation in (II.58) and the second part to 
the correction term. In (III.17) the linear term in @, Should be considered as 
a rather small correction and oY” is therefore determined only in the Gaussian 
approximation. The three parts are explicitly given below and our task is to 
put them into at form convenient for numerical calculations. 


Ag (Vn\"*# T' (n +8) [m _ 
uw’ | 


ue 


Eien) (n+l) 


of” 0) A, en@* Dana (2 Vn); (111.18) 


pow Al ey Teh ln, 
i (22) nT (n+1) mu’ 


3 
+ (30 3) x? + (nn —2 Demag (a Vn) — 
—aVn(na®+2n)D_nay(aVn);  (IIT.18’) 


Aot (Vn)"**T (n +38) [ See i Bee Bae 
gir ae nT el [7 7(0) Ae | é eet 


) T 9 Ay — — 2n 4n .— 
+ (Gat ght A8)| Ve D- cry (ela) Pelee yale 


+ (24242523) = (Se 2) pAs| Dens y le Vn)| PONTTLAR 4} 


The derivation of the formulae (III.18) is rather lengthy but elementary. The 
relation 


3 
(n 3) Pe Te Dnak @) dienyl) (111.19) 


between the Weber—Hermite functions has been used to express the cross sec- 
tion in only two of them and (m/M)y(0) Ay has been exchanged for (x— A,). 
Finally the terms are suitably rearranged. 


Approximate formula for o{'3 


In further simplifying the expression for oj the approximate expression for 
D_(n+y (2Vn) derived in Appendix III, is inserted into (III.18) and the pro- 


portionality factor depending on , which then appears, is approximated as 
follows: 
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Viel (nea) ee (III.20) 
| a a a 8n 


The terms are finally arranged in powers of 1 /n and we then get the following 


expression for o{"}: 


lm Ab +h) be ee | | 
te 2x3 (1 +327) z ry 1 (X) | 


a and z depend on the vibrations of the crystal according to 
m —_ 
v= (3 9(0)+8) Ao, 


m ne eo 1 a 2\" 


and h, (x) is defined by 


La[(l+ia?)?+ia2] 1 
a) 1+12° 8 


h, («) + @, (#). (111.23) 


y,() is given in Appendix III. The contribution by this term is only a few 
per cent and may in general be neglected. 

To get an idea of the validity of the approximations made in the derivation 
we will consider the two special cases x=0 and #= oo in more detail. In these 
eases (III.15) is easily integrated exactly and a comparison with the approxi- 
mate values obtained from (III.21) can be made. 

x=0: The exact and the approximate expressions differ only by the fol- 
lowing proportionality factor 


approx. values’ (n +1) ene Oe (IIT.24) 
exact value (/2n)"*#T n+3) 48 n . 


which is given numerically in Table 1 for some values of n. 


Table 1. 
n | il | 2 | 33 | 4 
ratio | 0.966 | 0.992 | 0.996 | 0.998 


a=cc: In this case the two formulae differ by the following factor 
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approx. value _ VnT (n+1) 1 3 (11.25) 
exact value T (n+) ’ 


and the numerical values are given in Table 2 


Table 2. 
esc sal geal 
ratio | 1.035 | 1.011 | 1.005 1.002 


Approximate formula for o\} 


The expression for o{ in (III.18’) has been approximated in essentially the 
same way as oj. The derivation is more complicated, however. To get a con- 
sistent approximation for all values of x the 1/n” term in the asymptotic ex- 
pansion of the Weber—Hermite functions must also be considered (see Appen- 
dix III). To simplify the final expression some terms that are small for all values 
of a have been neglected. The result may be written in the following form: 


AG (sri)! he 
2n* (1+1427)8 


1 | ; 
2" Vo jn (%) + Jp (@) + 79 (a) + | | (111.26) 


| oe = 
| 


The functions g,, g, and g, are given below. 


o-(Vne-9). 
ce 21+} aa) (+4 :) - 62 (Vise 147-3), 


4 )0Rs oe anl ee dee 
ial 323 — L+: +32%+1)- 
Js (@) 2 | ipa [ oe 2.) ( ) 
coe 259, =] 
(144) ae 
2 
_ Ae whe 5.1, 
16 32 128 
=3-27x7 +--+, «e>l1 


We notice that as x increases (8->0) and the high phonon terms contribute 
more and more to the cross section both g, and g, tend to zero. Owing to 


357 


A. SJOLANDER, Multi-phonon processes in neutron scattering by crystals 


this fact the Gaussian approximation will be more accurate for small values 
of 8. Such phonon terms as are not at all approximated by the Gaussian ap- 
proximation (very large values of x) do not contribute appreciably. 


Approximate formula for of” 


The derivation is made as for o{"j, so we only write the final result. 


n 1 
i A, 228 (+142) 2 enters ntorstmiors], (111.30) 


lo we = £1. te ter ~, (2 
potey= (5-5-4 4--) +p A,(G+5+5+~)-pbs (Ge), e< 


=(—84 2+ +--)+BA,(23a° 2+ ---)— pA? (—jet+--),emd 
p, and p, are given only for small and large values of 2. 


Summation over the phonon terms 


The summation over the phonon terms is now easily made. Including the 


one-phonon term we get the following expression for the inelastic cross section 
(cf. formulae (III.6) and (III.9)). 


TIT.31 
‘tetra [% 40,084] : ) 
0 


with 


_ Ad (144%)? -h2 Y - 1 
Sra ad (1+127)4 [1+ Foe) e + a +9 J3) In (=) de 


= z ; 
+ V2 Jy (1 =| , (117.31 ) 


= @ 1 Sr 
7 2x3 (1+32%) A, Eres rere In (=)| : (111.31 ) 
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For B—0 (I-—-) c, and c, tend to 


coo (st), fe Siew cole = ie 
On wat \2MB)? nts \omp) °° 8 a’ 
and (III.31) thus approaches the expression given in (II.17’). 


Numerical comparisons 


Accurate calculations of the total cross section have been made using Walker’s 
frequency distribution for Al.(!) The contribution from each phonon term is 
determined and these values are compared with those obtained from the ap- 
proximate expressions (III.21) and (III.26). Two cases are considered, one with 
T=300°K and one with 7’ =74.8°K (see Tables 3a and 3b). The separate 
terms in the m/M expansion are also calculated from formula (III.13’) and 
these values are in Table 4 compared with the values which are obtained from 
Placzek’s expression for a Debye crystal (see ref. [7]). The Debye temperature 
is then defined from the Debye-Waller factor. For (m/M)¥7(0)A, we have 
taken the value 0.5.(?) 


Table 3a. (T’'=300°K) 


n | ofP exact | of | ofr of 
1 0.05593 0.04849 0.05183 
2 0.02046 0.01950 0.02044 
3 0.00813 0.00784 0.00812 
4 0.00323 0.00315 0.00323 
5 0.001281 0.001267 0.001288 
6 Shas 0.00051 0.00051 
7 0.00020 0.00020 
8 ares 0.00008 0.00008 

| 0.0898 | 0.0810 | 0.0857 

Table 3b. (T'=74.8°K) 

0,502, for instance, means 0.00503. 
n oa” (exact) of | 03 af of? 
l 0,502 0.459 0,478 
2 03,230 0,224 03232 
3 04115 0,111 0,114 
4 0,575 0,563 0,579 
5 0,288 0,281 0,288 


(4) This frequency spectrum is adopted only to get a realistic one for our comparisons. 
(2?) This case should for 7 =300°K correspond to a hypothetical crystal with the mass 
ratio M/m=5.46 and with the same frequency spectrum as Al. 
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Table 4. (T'=300°K) 


de 0 we 
ue )| On 


n | exact | Placzek’s approx.(?) 
1 +0.1017 +0.1019 
2 — 0.0062 — 0.0071 
3 — 0.0069 — 
+ + 0.0012 — 
| 0.0898 | 0.0948 


(1) Placzek does not give more than two terms in the expansion. 
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APPENDIX I 


The scattering cross section for neutrons colliding with free nuclei with a 
prescribed momentum distribution (1) 


Consider first a neutron colliding with a free nucleus having the velocity u 
before the collision. As previously the neutron energies before and after the 
collision are denoted by iw, and hw respectively. The energy-momentum laws 
then give 

2 


syrt («i= 0 (A-I.1) 


W— Wot 


and the scattering cross section may be written 


do ack hx 3 
Th Edo 45a +x-a), (A-I.2) 


a being the bound scattering length for the nucleus. The Fourier-integral repre- 
sentation of the 6-function is introduced and the cross section is averaged over 
the prescribed velocity distribution of the nucleus, adopting the same distribu- 
tion as found for the nucleus in the crystal considered. We get 


cies, | an een hid or 
(sea0),7° hoy | ex i (« cn 5) ] <exp [it(x-u)]>;dt. (A-1.3) 


In the harmonic approximation for the crystal vibrations this velocity distri- 
bution is a Gaussian one and we can then write 


: t? 
<exp [it (%+ it)])p= exp <5 <0) | (A-L.4) 
By differentiating <u(t’) u(t)>7, explicitly given in (1.8), we can further show that 
<(x% +1)» ee St Poe eo (A-1.5) 
Tr NM = s Ss M 
and putting this into (A-I.3) the scattering cross section is written 
@o pai. hie | 2h 
ee A , = Ades ten teeeeest\ lat: 
(55-) rae | |i ( Oot oaF t} exp 27 (K+ &+x) 
(A-I.6) 


This expression may now be transformed into the form given in (II.7). 


(4) Se also: G. F. von Dardel, Phys. Rev., 94, 1274 (1954), part. D. 
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APPENDIX II ° 


Tables of y(0), Ay, Yr» V2 and related quantities, obtained in the 
Debye approximation 


= BLD, where wp=the Debye frequency. 
2k,T 

Bp | yp (0) | Ao,p | By yp (0) exp (—fi/2) V1 | Ve 
0.0 co 0.577 0.000 co 0.000 — 0.050 
0.1 30.030 0.577 0.052 29.998 

0.2 15.066 0.575 0.116 14.965 — 0.023 — 0.049 
0.3 10.099 0.571 0.173 9.949 

0.4 7.632 0.566 0.231 7.431 — 0.046 — 0.045 
0.5 6.165 0.560 0.289 5.912 

0.6 5.197 0.554 0.348 4.892 — 0.067 — 0.040 
0.7 4.515 0.546 0.406 4,157 

0.8 4.010 0.537 0.464 3.600 — 0.088 — 0.032 
0.9 3.624 0.528 0.523 3.161 

1.0 3.321 0.518 0.582 2.804 — 0.107 — 0.024 
Lak 3.078 0.508 0.640 2.508 

1.2 2.879 0.497 0.699 2.255 — 0.123 — 0.014 
1.3 2.715 0.486 0.757 2.038 

1.4 2.577 0.476 0.815 1.848 — 0.138 — 0.004 
1.5 2.461 0.465 0.873 1.681 

1.6 2.362 0.455 0.931 1.532 — 0.150 + 0.005 
Vai 2.276 0.445 0.987 1.398 

1.8 2.203 0.435 1.044 1.278 — 0.160 +0.014 
1:9 2.138 0.425 1.099 1.169 

2.0 2.082 0.416 1.154 1.070 — 0.167 + 0.021 
2.5 1.885 0.376 1.411 0.697 

3.0 1.771 0.345 1.635 0.466 — 0.175 + 0.036 
3.5 1.701 0.325 1.823 0.323 

4.0 1.654 0.306 1.978 0.234 — 0.162 + 0.028 
4.5 1.622 0.293 2.104 0.177 

5.0 1.599 0.283 2.208 0.140 — 0.147 +0.017 
5.5 1.582 0.276 2.292 0.114 

6.0 1.569 0.270 2.361 0.097 — 0.136 + 0.006 
oo 1.500 0.236 2.825 0.027 — 0.094 — 0.025 
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Tables of (0), Ay and V5, obtained in the Debye approximation. 


(ooo 


Bp 


SD et et ee best ea oes Feel ety Feed 
SMe MmHIAnrwWNeE © 


YD (0) 


Ao,D 


0.577 
0.577 
0.577 
0.576 
0.575 
0.574 
0.573 
0.571 
0.569 
0.567 
0.565 
0.562 
0.560 
0.557 
0.554 
0.550 
0.547 
0.543 
0.540 
0.536 
0.532 
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Ya 
— 0.0500 
— 0.0498 
— 0.0497 
— 0.0493 
— 0.0488 
— 0.0481 
— 0.0473 
— 0.0463 
— 0.0453 
— 0.0440 


— 0.0427 
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Tabb 
© 2” exp (— 27/2m 
F(z, x)= — an 
neg tl (227) 
x= 0.0, | 
be 
a 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 | 


0.0 | 0,145 0,596 0,138 0,253 0,406 0,603 0,846 0.114 0.149 0.190 08 
0.1 | 0,145 0,595 0,138 0,252 0,405 0,601 0,844 0.114 0.148 0.189 O¢ 
0.2 | 0,144 0,591 0,137 0,250 0,402 0,597 0,838 0.113 0.147 0.188 0.8 
0.3 | 0,142 0,583 0,135 0,247 0,398 0,590 0,828 0.112 0.146 0.186 0. 
0.4 | 0,139 0,573 0,133 0,243 0,391 0,580 0,815 0.110 0.143 0.183 O£ 
0.5 | 0,136 0,561 0,130 0,238 0,383 0,568 0,798 0.108 0.140 0.179 02 
0.6 | 0,133 0,546 0,127 0,232 0,373 0,554 0,777 0.105 0.1387 0.175 8 
0.7 | 0.128 09529 0,123 0,225 0,362 0,537 0,754 0.102 0.133 0.170 0.2 
0.8 | 05124 0,510 0,118 0,217 0,349 0,518 0,728 0,983 0.129 0.164 0.2 
0.9 | 0,119 0.489 0,113 0,208 0,335 0,498 0,700 0,945 0.124 0.158 0.1 
0,113 0,467 0,108 0,199 0,320 0,476 0,670 0,905 0.118 0.151 0,1 
0,107 0,443 0,103 0,189 0,305 0,453 0,638 0,862 0.113 0.144 1 
0,102 0419 0,973 0,179 0,289 0,430 0,605 0,818 0.107 0.137 0.1 
0,954 0,394 0,916 0,168 0,272 0,405 0,571 0,772 0.101 0.130 0.1 
05892 0,369 0,858 0,158 0,255 0,380 0,536 0,726 0,953 0,122 O.1 
05831 0,344 0,800 0,147 0,238 0,356 0,502 0,680 0,893 0.114 0.1 
03769 0,319 0,742 0,137 0,221 0,331 0,467 0,633 0,833 0.107 OJ 
03709 05294 0,685 0,126 0,205 0,306 0,433 0,588 0,773 0,993 0.1 
0;650 0,270 0,630 0,116 0,189 0,282 0,400 0,543 0,715 0,919 0.1 
03594 0,247 0,576 0,107 0,173 0,259 0,367 0,499 0,659 0,848 0.1 
03539 0,224 0,525 0,970 0,158 0,237 0,336 0,457 0,604 0,778 0,9 
03487 0,203 02475 0,880 0,143 0,215 0,306 0,417 0,552 0,712 0,9 
09438 02183 02429 0,795 0,130 0,195 0,278 0,379 0,502 0,648 0,8 
03392 0,164 0,385 02715 0,117 0,176 0,251 0,343 0,455 0,588 0,7 
05349 0,146 02344 02640 0,105 += 0,158 04226 ~=— 0,309 0,410 0,532 0,6 
0,310 02130 0,306 02570 0,935 0,141 0,202 0,277 0,369 0,479 0,¢ 
03273 0,115 0,271 02506 += 02831 0,126 += 04180 04248 0,331 0,429 045 
03240 0,101 0,239 03447 0,736 0,112 0,160 0,221 0,295 0,384 044 
05209 03883 0,210 02393 02648 0.986 0,142 0,196 0,262 0,342 0,4 
03182 0,770 0,183 02344 0,569 0,868 0,125 0,173 0,232 0,304 043 
0,158 0,668 02159 0,300 02498 0,760 0,110 04,152 0,205 0,269 0,3 
0,136 05577 0,138 0,261 02433 0.644 ~— 02962 0,134 0,180 0,237 048 
0,116 0,496 0,119 0,226 0,376 ~— 02578 02839 0,117 0,158 0,208 042 
04992 03424 0,102 0,194 0,325 02501 02729 0,102 04138 0,183 092 
04842 0,362 0,874 0,167 0,280 0,432 0,632 0,886 0,120 0,160 0,2 
0,711 03307 03744 +0143 0,240 04372 0546 = 04768 0,105 0,189 04d 
04598 05259 03631 0,121 0,205 04319 02470 0,663 0,907 04,121 4] 
0,501 03218 03533 0,103 04175 04273 0g403 04571 0,784 0,105 04] 
-8 | 04418 05183 05449 0,871 02148 0,233 034504491 04676 += 0,907 04 
3.9 | 04346 05152 04376 0373403126 0g198 04295 0,421 04581 0478304 

04286 03127 05314 03616 0,106 0,168 04251 0,360 05499 04675 Oe8 
04235 03105 04262 03516 0589204142 0.213 04307 0427 04580 097 
04193 04864 05217 0,431 05748 05120 04181 0,261 04365 0,498 Og 
4.3 | 04157 04709 05180 0,858 0,626 0.101 04153 0,222 04,312 0426 Oat 
4.4 | 04128 0,581 05148 05297 03523 0,845 0,129 0,188 04,265 04365 Od 
4.5 | 04103 04473 05122 05246 05436 05708-04109 04159 05226 04311 Od 


WH SPP PNNNNNN Epo pee poe 


URDU PR WN HOOBDUAMURWNHHOOHDUDMARWHHS 


g2 92 G2 G2 O Oo 


oP 
Sreo 


4.6 0,592 05913 0,135 05192 0,266 043 
4.7 05494 05766 0114 0,163 0,226 045 
4.8 0,412 0,642 0,957 0,138 0,192 043 
4.9 0,343 05538 0,805 05116 0,164 048 
os 0,285 0,450 0,677 0,984 04139 0g] 
5.2 
5.3 
5.4 
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91 0.353 0.424 0.504 0.593 0.694 0.806 0.932 1.071 1.226 1.398 
0.353 0.423 0.502 0.592 0.692 0.805 0.930 1.069 1.224 1.396 
0.350 0.420 0.500 0.588 0.688 0.800 0.924 1.063 1.221 1.387 
85 0.346 0.416 0.494 0.582 0.681 0.791 0.915 1.052 1.205 
0.341 0.409 0.487 0.574 0.671 0.780 0.902 1.037 1.188 
76 0.335 0.402 0.477 0.563 0.659 0.766 0.886 1.019 1.167 
69 0.327 0.392 0.466 0.550 0.644 0.749 0.866 0.997 1.142 
tL 0.318 0.381 0.454 0.5385 0.627 0.729 0.844 0.971 1.113 
53 0.307 0.369 0.489 0.519 0.607 0.707 0.818 0.943 1.081 
44 0.296 0.356 0.424 0.500 0.586 0.683 0.791 0.911 1.046 
34 0.284 0.342 0.407 0.481 0.564 0.657 0.761 0.878 1.007 
23 0.272 0.327 0.390 0.460 0.540 0.630 0.730 0.842 0.967 1,106 
12 0.258 0.311 0.371 0.439 0.515 0.601 0.697 0.805 0.925 1.059 
01 0.245 0.295 0.352 0.417 0.490 0.572 0.663 0.766 0.881 1.009 
0.231 0.279 0.333 0.394 0.463 0.541 0.629 0.727 0.837 0.959 
0.217 0.262 0.313 0.371 0.437 0.511 0.594 0.687 0.791 0.908 
0.203 0.246 0.294 . 0.349 0.411 0.480 0.559 0.647 0.746 0.857 
0.190 0.229 0.274 0.326 0.384 0.450 0.524 0,607 0.701 0.805 
0.176 0.213 0.255 0.304 0.358 0.420 0.490 0.568 0.656 0.755 
0.163 0.197 0.237 0,282 0.333 0.391 0.456 0.530 0.612 0.705 
0.150 0.182 0.219 0.261 0.308 0.3862 0.423 0.492 0.569 0.656 
0.138 9.167 0.201 0.240 0.284 0.335 0.391 0.456 0.528 0.609 
0.126 0.153 0.185 0.220 0.261 0.308 0.361 0.421 0.488 0.564 
0.115 0.140 0.169 0.202 0.240 0.283 0.332 0.387 0.450 0,520 
0.104 0.127 0.154 0.184 0.219 0.259 0.304 0.355 0.413 0.479 
0,943 0.115 0.140 0.167 40.199 0.286 0.278 0.325 0.379 0.439 
0,851 0.104 0.126 0.152 0.181 0.215 0.258 0.296 0.346 0.402 
0,765 0,988 0.114 0.137 0.164 0.195 0.230 0.270 0.315 0.367 
0,685 0,842 0.102 0.124 0.148 0.176 0.208 0.245 0.287 0.334 
0,612 0,754 0,919 0.111 0.183 0.159 0.188 0.222 0.260 0.303 
437 0,545 0,673 0,822 0,995 0.120 0.143 0.169 0.200 0.235 0,275 
387 0,484 0,599 0,733 0,889 0.107 0.128 0.152 0.180 0.212 0.249 
342 0,429 0,532 0,652 0,793 0,957 0.115 0.137 0.162 90.191 0.225 
301 0,379 0,471 0,579 0,705 0,853 0.103 0.122 0.145 0.172 0.202 
265 0,334 0,416 0,512 0,626 0,759 0,914 0.109 0.130 0.154 0,182 
232 0,293 0,366 0,453 0,555 0,674 0,813 0,975 0.116 0.138 0.163 
203 0,257 0,322 0,399 0,490 0,597 0,722 0,868 0.104 0.124 0.146 
1 ae 0,225 0,283 0,351 0,432 0,528 0,640 0,772 0,925 0.110 0.131 
154 0,197 0,248 0,308 0,381 0,466 0,567 0,685 0,822 0,982 0.117 
134 0,171 0,216 0,270 0,335 0,411 0,501 0,607 0,730 0,874 0.104 
116 0,149 0,189 0,237 0,294 0,362 0,442 0,537 0,647 0,777 (0,928 
101 0,130 0,165 0,207 0,258 0,318 0,390 0,474 0,573 0,690 0,826 
870 0,112 0,143 0,181 0,225 0,279 0,343 0,418 0,507 0,612 0,734 
751 0,973 0,124 0,157 0,197 0,245 0,301 0,369 0,448 0,542 0,651 
1647 0,841 0,108 0,137 0,172 0,214 0,265 0,824 0,395 0,479 0,577 
557 0,727 0,936 0,119 0,150 0,187 0,232 0,285 = 0,349 0,423 0,511 
A479 «0,627 = 0,810 = 0,103 04131 0,164 90,2038 0,251 0,807 0,874 0,453 
411 0,540 0,700 0,897 0,114 0,143 0,178 0,220 0,270 0,330 0,400 
353 0,465 0,605 02778 0,989 0,125 0,156 0,193 04237 0,290 0,353 
302 0,400 0,522 0,673 0,859 0,109 0,136 0,169 0,208 0,256 0,312 
259 0,344 0,450 0,583 02745 0,944 0,119 0,148 0,183 0,225 0,275 
0,821 0,104 0,129 0,160 0,198 0,242 
0.714 0,902 0,113 0,141 0,174 0,218 
0,620 0,785 0,986 0,123 0,152 0,188 
0,538 0,683 0,861 0,108 0,134 0,165 
ne 
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(contd.) 


2\2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 


0.0 | 1.589 1.799 2.032 2.288 2.571 2.883 3.226 3.604 4.019 4.475 
0.2 | 1.577 1.786 2.017 2.272 2.553 2.862 3.203 3.579 3.991 4.445 
0.4 | 1.540 1.745 1.972 2.222 2.497 2.801 3.136 3.504 3.910 4.355 
0.6 | 1.482 1.680 1.899 2.141 2.408 2.703 3.027 3.384 3.778 4.210 
0.8 | 1.405 1.594 1.803 2.034 2.289 2.571 2.882 3.224 3.600 4.017 
1.0 |~-1.312 1.490 1.687 1.905 2.146 2.412 2.706 3.031 3.388 3.783 
2 | 1.207 1.372 1.556 1.759 1.984 2.233 2.508 2.812 3.147 3.517 
4 | 1.095 1.247 1.415 1.602 1.810 2.040 2.294 2.575 2.887 3.231 
6 | 0.080 1.117 1.270° 1.441 1.630% "1.840 2.072" 2.330" 2.616 = se gae 
8 | 0.865 0.988 1.126 1.279 1.450 1.639 1.850 2.084 2.344 2.631 
0 | 0.754 0.863 0.985 1.122 1.274 1.4444 1.633 1.843 2.077 2.336 
2 | 0.649 0.745 0.853 0.973 1.107 1.258 1.426 1.613 1.821 2.053 
2.4 | 0.553 0.636 0.730 0.835 0.952 1.084 1.232 1.397 1.581 1.787 
2.6 | 0.465 0.537 0.618 0.709 0.811 0.926 1.055 1.199 1.361 1.542 
2.8 | 0.388 0.449 0.518 0.596 0.685 0.784 0.895 1.021 1.162 1.320 
3.0 | 0.321 0.372 0.481 0.498 0.573 0.658 0.754 0.862 0.984 1.121 IL: 
3.2 | 0.263 0.306 0.356 0.412 0.476 0.548 0.630 0.723 0.827 0.945 I 
3.4 | 0.214 0.250 0.291 0.339 0.392 0.454 0.523 0.602 0.691 0.792 0! 
3.6 | 0.172 0.202 0.237 0.277 0.322 0.373 0.432 0.498 0.574 0.660 0,’ 
3.8 | 0.138 0.163 0.192 0.225 0.262 0.305 0.354 0.411 0.474 0.547 Ok 
4.0 | 0.110 0.131 0.154 0.181 0.213 0.248 0.289 0.336 0.390 0.451 0. 
4.2 | 0,876 0.104 0.124 0.146 0.171 0.201 0.235 0.274 0.319 0.371 © 
4.4 | 0,693 0,827 0,984 0.117 0.138 0.162 0.190 0.223 0.260 0.303 0: 
4.6 | 0,545 0,654 0,781 0,930 0.110 0.130 0.154 0.180 0.211 0.247 0: 
4.8 | 0,428 0,515 0,618 0,738 0,879 0.104 0.128 0.145 0.171 0.200 OG: 
5.0 | 0,334 0,404 0,487 0,584 0,698 0,831 0,986 0.117 0.138 0.162 0. 
5.2 | 0,260 0,316 0,383 0,461 0,553 0,661 0,787 0,934 0.111 0.131 O11 
5.4 | 0,202 0,247 0,300 0,362 0,436 0,523 0,626 0,745 0,885 0.105 11 
5.6 | 0,157 0,192 0,234 0,284 0,343 0,413 0,496 0,593 0,706 0,839 Og 
5.8 | 0,121 0,149 0,182 0,222 0,269 0,326 0,392 0,470 0,562 0,670 0, 
6.0 | 0,930 0,115 0,141 0,173 0,211 0,256 0,309 0,372 0,446 0,533 0, 
6.2 | 0,714 0,887 0,109 0,134 0,164 0,200 0,243 0,293 0,353 0,423 04! 
6.4 | 0,547 0,682 0,845 0,104 0,128 0,156 0,190 0,230 0,278 0,335 04 
6.6 | 0418 0,523 0,651 0,805 0,992 0,122 0,149 0,181 0,219 0,264 0, 
6.8 | 0,318 0,400 0,500 0,621 0,768 0,945 0,116 0,141 0,172 0,208 04: 
7.0 | 0,242 0,305 0,383 0,478 0,593 0,732 0,900 0,110 0,134 0,163 041 
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3.7 


6.131 6.795 7.523 
6.092 6.752 7.476 
5.974 6.624 7.336 
5.785 6.416 7.110 
5.530 6.139 6.807 
5.222 5.802 6.439 
4.872 5.418 6.020 
4.492 5.003 5.566 
4.096 4.568 5.089 
3.695 4.127 4.606 
3.299 3.692 4.128 
2017 © 3.272. « 3.665 
2.557 2.873 - 3.226 
2.222 2.503 2.816 
1916 2.164 2.441 
1.641 1.857 2.100 
1395 1.584 1.796 
Hise «41.848 «= «1.527 
0.992 1.132 1.291 
0.830 0.950 1.086 
0.691 0.793 0.909 
0.572 0.659 0.758 
0.472 0.546 0.629 
0.388 0.450 0.520 
0.318 0.370 0.429 
0.260 0.303 0.352 
0.211 0.247 0.288 
0.171 0.201 0.235 
0.138 0.163 0.191 
Gili, 0.131 0,155 
0,895 0.106 0.125 
0,717. 0,850 0.101 
0,573 0,681 0,808 
0,456 0,544 0,647 
0,362 0,433 0,517 
0,287 0,344 0,412 
bey. 0,973 . 0,328 
0,179 0,216 0,260 
0,141 0,170 0,205 
0,111 0,134 0,162 
0,865 0,105 0,128 


3.9 


9.199 
9.143 
8.977 
8.709 
8.349 
7.912 
7.413 
6.870 
6.300 
5.720 
5.144 
4.585 
4.053 
3.554 
3.094 
2.676 
2.300 
1.965 
1.670 
1.413 
1.189 
0.997 
0.832 
0.692 
0.573 
0.473 
0.389 
0.319 
0.261 
0,212 
0.173 
0.140 
0.113 
0,909 
0,730 
0,585 
0,468 
0,373 
0,296 
0,235 
0,186 
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4.0 4.1 4.2 
10.16 41:29 12.37 
103 OS HIET SINS 12730 

9.920 10.95 12.09 

9.628 10.64 11.74 

9.236 10.21 11.28 

8.760 9.692 10.72 

8.216 9.100 10.07 

7.623 8.454 9.368 
7.001 7.774 8.627 
6.366 7.080 7.869 
5.736 6.390 7.113 
5.122 5.717 6.376 
4.536 5.073 5.669 
3.987 4.468 5.003 
3.479 3.908 4.385 
3.016 3.395 3.818 
2.598 2.932 3.305 
2.226 2.518 2.845 
1.896 2.151 2.436 
1.608 1.828 2.076 
PAT eailbi ten aa 1e! 
M141) (418035,0 1,488 
0.955 1.094 1.251 
0.796 0.914 1.049 
0.661 0.761 0.875 
0.547 0.632 0.728 
0.451 0.522 0.604 
0.371 0.430 0.499 
0.304 0.354 0.411 
0.248 0.290 0.337 
0.202 0.237 0.276 
0.164 0.193 0.225 
0.133 0.156 0.183 
0.107 0.127 0.149 
0,865 0.102 0.120 
0,694 0,823 0,972 
0,556 0,661 0,783 
0,445 0,530 0,629 
0,355 0,423 0,504 
0,282 0,338 0,403 
0,224 0,269 0,321 


13.64 
13.56 
13.33 
12.96 
12.46 
11.84 
11.14 
10.37 
9.567 
8.740 
7.913 
7.105 
6.330 
5.598 
4.917 
4.291 
3.723 
3.212 
2.757 
2.355 
2.003 
1.696 
1.430 
1.201 
1.005 
0.838 
0.697 
0.577 
0.476 
0.392 
0.322 
0.263 
0.215 
0.175 
0.142 
0.115 
0,926 
0,746 
0,599 
0,480 
0,384 
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APPENDIX III 
Approximate expressions of the Weber—Hermite functions (?) 


The Weber-Hermite functions D_(n+3) (x Vn) satisfy the following differential 
equation: 
2 


y"’ (a) —n® (1 + =) y (x) =0, (A-IIT.1) 


and as we are especially concerned with rather large values of n(m>2) we can 
use the WKB method to obtain approximate expressions for these functions. 
According to this method a new function w(x) is defined by 


y (x) = exp (n/ u (t) at) (A-ITT.2) 
0 


and this satisfies the equation below 


1 2 
~w=1 ae uw, (A-III.3) 


If w(x) is now expanded into powers of 1/n 


u(x) = 2 n © Uy (2). (A-ITT.4) 


u, are determined from the following recursion formulae: 


ug=1+32'%, 
Uo +2 uy u,=9, 
Uy +2 Uy Uy + uj=0, (A-IIT.5) 


Uz +2 Uy Ug + 2 U, Uy =0, 


and it is a straightforward matter to derive the solutions. The first few of 
these are given below. 


(*) The properties of these functions are briefly discussed in e.g. Erdelyi, Magnus, Ober- 
hettinger, and Tricomi, Higher Transcendental Functions, vol. 2, pp. 116-130, New York 
1953. See also: N. Schwid, Trans. Amer. Math. Soc., 37, 339 (1935). 
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4 
d 2\—4 
Nae In 144) ; 
eee! elas 2\—3 A-IIL.6 
= F 35 [3 (144) -5(1+4) |: 


The original function y(#) is then obtained in the form 


const 3 £2\4 1 ; 1 : 
0 0 0 
(A-III.7) 


The + signs correspond to two independent solutions, and we must specify 


both the sign and the constant factor to get an approximation for D_(n+4) (x Vn). 
For this purpose we let x->co, where (A-III.7) tends to 


- 2 
y («)-> const V2 e° (Vn)!*" (x Vn)-#*" exp (+72) (A-IIT.8) 
. i ) 1 C) 
with c= Bat fagydie ts | u( drs ves 
0 0 


This expression should be compared with the corresponding one for D_¢n +4, (a Vn) (2) 
D_ n+) (a Vn)—>(aVn)-@t® e284 (A-ITI.9) 


leading us to conclude that the lower sign should be chosen and that the 
constant factor must be 


const = 2>* e~¢ (Vn)- +», (A-IIT.10) 


To obtain a more convenient expression we finally expand in powers of 1/n, 
where we get 


_ acme [=n(2-+/0)] +20 e+ 


(A-IIT.11) 


(1) See: Erdelyi et al., loc. cit., p. 122. 
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l a?\t gl]? a\t 2 
where j= -5|(1+4) -{| —In (1+) ms , 


(A-IIL.11’) 


A corresponding approximate expression for D_ n+ p (x Vn) is obtained from 
(A-III.11) by using the relation 


al d 1 
(» + 3) D_(n+3 (2) = — exp is i?) mF jexp (; 2) D_wn+)) (. (A-ITT.12) 


We get 


ee: 2a 


with ty (2) = (2) +5 aA’ (A-IIT.13’) 
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